Accepted by the Astrophysical Journal 

Preprint typeset using I^FfjjK style cmulatcapj v. 2/16/10 



ON THE REDSHIFT-EVOLUTION OF THE LYMAN- ALPHA ESCAPE FRACTION AND THE DUST 

CONTENT OF GALAXIES 

Matthew Hayes 1 , Daniel Schaerer 1,2 , Goran Ostlin 3 , J. Miguel Mas-Hesse 4 , Hakim Atek 5 , Daniel Kunth 6 

Accepted by the Astrophysical Journal 

ABSTRACT 

The Lya emission line has been proven a powerful tool by which to study evolving galaxies at the 
highest redshifts. However, in order to use Lya as a physical probe of galaxies, it becomes vital to 
know the Lya escape fraction (/ e ^ a ). Unfortunately, due to the resonant nature of Lya, may 
vary unpredictably and requires empirical measurement. Here we compile Lya luminosity functions 
between redshift z=0 and 8 and, combined with Ha and ultraviolet data, assess how evolves 
with redshift. We find a strong upwards evolution in f^. a over the range z — 0.3 — 6, which is well-fit 
by the power-law f^/ c a (x (1 + z)« with £ = (2.57lto'.i2)- This predicts that / L s y Q should reach unity at 
z = 11.1. By comparing f^. a and Eb-v in individual galaxies we derive an empirical relationship 
between f^ a and Eb-v, which includes resonance scattering and can explain the redshift evolution 
of f^/ c a between z — and 6 purely as a function of the evolution in the dust content of galaxies. 
Beyond z ~ 6.5, f^. a drops more substantially; an effect attributed to either ionizing photon leakage, 
or an increase in the neutral gas fraction of the intergalactic medium. While distinguishing between 
those two scenarios may be extremely challenging, by framing the problem this way we remove the 
uncertainty of the halo mass from Lya-based tests of reionization. We finally derive a new method 
by which to estimate the dust content of galaxies based purely upon the observed Lya and UV LFs. 
These data are characterized by an exponential with an e-folding redshift of ps 3.5. 
Subject headings: Galaxies: evolution — Galaxies: high-redshift — Galaxies: luminosity function, 
mass function — Galaxies: star formation — dark ages, reionization, first stars 



1. INTRODUCTION 

Surveys targeting the Lyman-alpha emission line (Lya) 
show unique profitability for examining the formation 
and evolution of the galaxy population between redshift 
z « 2 and > 7. Lya has been exploited by many teams 
and the combined catalogues w ould currently include 
over two thousand entries fe.g. [Venemans et al.l l2002t 
Hu et all 120041 Ivan Breukelen et al.1 120051: iWang et al.1 
20051: IShimasaku et al.l I2006t IGronwall et al.l 1 20071: 
Ouchi et al.ll2008l : lNilsson et al.ll200a iGuaita et al.ll20ld 



Cassata et alj|2010t lHaves et al.ll2010aD . Wherever such 
large samples are available, the temptation is strong to 
use their statistical power to examine as many physical 
properties of the galaxy population as possible. This, 
however, requires that the numbers one has at hand are 
in some way a physical reflection of those underlying 
properties; to first order the luminosity (and/or equiv- 
alent width for emission lines) must be related to its 
intrinsic value. For surveys that target the restframe 
ultraviolet continuum (UV) this is simply a matter of ap- 
plying a dust correction. However, the resonant nature 
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of the Lya line means its radiat ion transport becomes an 
invol v ed and detailed p r oblem (|Osterbroc k 1962; Adams! 
19721 iHarringtod 11^731: iNeufeldl 119901 I Arm et alj|2003t 
Verhamme et al.1 120061 : iTasitsiomil 120061 : lLaursen et al.1 
2009( 1. This further implies that the escaping fraction of 



photons (/ e L s f ) may not be assumed, is liable to evolve 
strongly with an evolving galaxy population, and must 
be measured empirically. Pursuing this line of inquiry, 
the evolution of can therefore provide us with inde- 
pendent estimates of how various properties of the galaxy 
population evolve over cosmic time. 

Since Lya photons scatter in neutral hydrogen (Hi) 
until they either escape or are absorbed by dust grains, 
most fundamentally the radiation transport depends 
upon the Hi content, its geometry and kinematics, and 
the dust content and distribution. Regrettably, with cur- 
rent observational facilities, the only one of these quanti- 
ties that can easily be estimated for large samples of high- 
redshift galaxies is the dust attenuation, which is typi- 
cally derived from the stellar continuum. Consequently 
the amalgamated effects of the remaining quantities, and 
how they affect f^/ c a , can only be assessed on a statistical 
basis. 

Lya surveys have been fruitful over the last decade, but 
it is only very recently that robust measurements 
have been made on st atistically meani n gful samples 



(Tverhamme et al.1 
120101: lHaves et al.l 



on statistically meaningful samples 
2001 I Atek etafl l200l iKornei et all 



2010al ). However at the current junc- 



ture, all of these studies estimate f^. a by different meth- 
ods, and are derived among samples compiled at various 
redshifts and filtered through differing selection func- 
tions. Thus synthesis of the results remains somewhat 
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difficult. Furthermore, there is no self-consistent study 
in the current literature of how f^. a evolves with redshift 
and it is this point that we take the first steps towards 
rectifying with the current article. We begin by com- 
piling various Lya, Ha, and UV datasets in § [21 which 
we use to estimate the redshift evolution of / c L s ^ a . We 
discuss the general trends and draw comparisons with 
other observational and theoretical methods in § [3] In 
§ |4] we investigate the effect of the one quantity that is 
relatively easy to measure - the dust content - and dis- 
cuss how it affects f^. a . In §[5] we discuss the trends with 
redshift in more detail and synthesize information from 
§ |31 and § [3] in order to make more detailed inferences 
about the evolution of the properties of the interstel- 
lar medium (ISM) of galaxies, the intergalactic medium 
(IGM), and the overall dust content. In §[5] we present 
a final summary. All data are scaled to a cosmology of 
(H ,ttM,aO = (70 km s^ 1 Mpc _1 ,0.3,0.7). 



2. METHOD: THE Lya ESCAPE FRACTION 
MEASUREMENTS 

2.1. Escape fraction calculations 

We now proceed to compile various estimates of f^/ c a 
as a function of redshift, but first we p r esent t he formal- 
ism. We continue with the lHaves et al.l ([2pi0a) definition 
of : the sample-averaged, "volumetric" escape frac- 
tion. This quantity is defined as the ratio of observed to 
intrinsic Lya luminosity densities (ph), derived by inte- 
gration over luminosity functions (LF), as in Equation [TJ 
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L-dL 



dL 



(1) 



where $(£) are the standard luminosity function^. Thus 
fe^. a is not simply a re-scaling of the LF by L (constantly 
scaling the es cape fraction of all ga laxies) or by $ (the 
duty cycle; see[Nagaminc et al. 2008 for examples of both 
of these methods). Instead, since f^. a is simply defined 
as the ratio of luminosity densities, it can be thought 
of as the fraction of Lya photons that escape from the 
survey volume, regardless of whether all galaxies show 
low f^. a , or whether only a fraction of galaxies are in 
the Lya emitting phase with high (see arguments 

in iTilvi et ani20CH) . By definition f^/ c a also includes 
any possible effect that the IGM may have on the Lya 
emission from galaxies. However, it is clear that the bulk 
of the evolution of f^. a with redshift found in this article 
can clearly not be attributed to variations of the IGM 
transmission. 

Where possible (i.e. z < 2.3) we make a direct com- 
parison between Lya and Ha. We apply the most appro- 
priate dust correction to Ha and mu ltiply by the case B 
recombination ratio of Lya/Ha=8.7 (Brocklc hurstill97l[l 
in order to obtain the intrinsic Lya. I.e. 
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7 LFs are typically parameterized by the Schcchtcr ( 1976) func- 
tion: *(L) • dL = ■ (L/L+) 01 ■ exp(L/L*) ■ dL/L+. 



where Eb-v must be the dust attenuation computed for 
the Ha emitting sample, and ^6563 the extinction coeffi- 
cient at the wavelength of Ha. Superscripts Int and Obs 
refer to the intrinsic and observed quantities. 

At z > 2.3 we are unable to obtain Ha LFs in or- 
der to use line ratios to estimate f^. a and instead the 
estimate is derived from the UV continuum. This is a 
less elegant method since the conversion between UV 
and Lya requires the assumption of a metallicity, initial 
mass function (IMF), and evolutionary stage. However, 
in light of the fact that higher-redshift Ha studies will 
remain impossible until the arrival of the James Webb 
Space Telescope, this is the only way to proceed. It is 
fortunate that there is no evidence that IMFs should dif- 
fer between Lya- and UV selected populations, although 
met allicities have been shown to be around 0.2 dex lower 
(e.g. lCowie et al.ll2010l) which translates into a difference 
of < 20 % in the intrinsic Lya/UV ratio (|Leitherer et al.l 
1999). For "normal" metallicities and IMFs, and assum- 
ing that on average star-formation is ongoing at equi- 
librium, this method is the same as taking the ratio of 
Lya/UV star- formation rate densities (/»*): 



> 2.3) - 
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where now Eb-v must be the extinction seen by the UV- 
selected population and fcrjv ^ s ^ ne extinction coefficient 
in the UV. 

The UV is of course not the only wavelength we can 
use for this experiment, but we choose to work exclu- 
sively with UV LFs since they (a) are so abundant in the 
literature, (b) have reasonably well-understood selection 
functions, and (c) span an appropriately large range in 
redshift. We adopt UV measurements at redshifts most 
appropriate to our compiled Lya data and dust attenua- 
tions derived from these samples t hemselves. We furthe r 
adopt the dust attenuation la w oflCalzetti et"aL1 (2000), 
and the SFR calibrations of iKennicuttl (|1998| ). These 
calibrations assume a stabilized star formation episode 
at a constant rate for longer than around 100 Myr, with 
a Salpeter Initial Mass function (mass limits between 0.1 
and 100 Mq), and a complete ionization efficiency (no 
leaking and no destruction of ionizing photons by dust). 
In general we assume that 'UV refers to the restframc 
wavelength of 1500 A, where the extinction coe f ficient 
computed from the relationship of .Calzctti et al. (2000) 
is 10.3. We want to emphasize that the definition of 
/esc a we are using for high redshift galaxies includes any 
effect that would decrease the number of observed Lya 
photons with respect to the number expected from the 
star formation rate derived from the UV continuum level. 
The leaking of ionizing photons, as we will discuss later, 
would therefore imply an f^ a value below unity, even 
if 100% of the Lya photons effectively produced in the 
galaxy are able to escape without being affected by res- 
onant trapping or destruction by dust. 

2.2. Limits of integration 

The goal of this study is to determine the total, vol- 
umetric escape fraction of a given volume, and ideally 
would include the very faintest systems. In practice this 
would require integration of the LFs down to zero, which 
depending on the observational limits of a given survey 
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and the redshift-dependent values of both L+ and a, may 
include large extrapolations (or may even be divergent). 
It is vital therefore, that our study employs lower in- 
tegration limits that are: (a) self-consistent between the 
populations; (b) include a sufficiently meaningful fraction 
of pl, and (c) are not dominated by over-extrapolation 
and uncertainties in the faint-end slope. 

At z — 2, 3 and > 4, several studies of th e pl.uv have 
been p ublis hed, and here w e adopt those of iReddv et al.l 
( 2008) and IBouwens et al.l (|2009aD . respectively. Both 
perform integrations down to 0.04LJ^j V and integrate to 
the same numerical lower limit at all redshifts. The lower 
limit is, of course, somewhat arbitrary but is designed to 
find a reasonable medium between including a large frac- 
tion of the total luminosity/ SFR density, and preventing 
(possible over-) extrapolation by integrating to zero. In 
this sense, it reflects the observational limits of the UV 
surveys. 

Admitting that this number is somewhat arbitrary, we 
adopt the same approach and use 0.04LJ =3 — oo as the 
range for all of integrations of the UV LF. For M^ =3 = 
—21.0 (AB), the corresponding lower luminosity limit is 
4.36 x 10 27 erg s" 1 Hz" 1 (unobscured SFR=O.6M vr -1 ). 
By adopting this limit, our results can easily be cross- 
checked agains t the available litera ture. At redshift 3 for 
the UV LF of Re ddv et al.l ()2008h . this range incorpo- 
rates 70% of an infinite integration under the LF. 

Deciding upon a lower limit for the Ha LF is more 
tricky, since it is difficult to know if we are extract- 
ing comparable samples of galaxies. There is no avail- 
able z = 3 Ha LF, but if we adopt that compiled at 
z = 2.2 in lHaves et al.l (|2010bl ). and set the lower limit 
to 0.04L*, we obtain a luminosity of 4.6 x 10 41 erg s _1 . 
This corresponds to much higher unobscured SFR than 
the lower UV limit at 3.5A/ Q yr -1 . However, the UV 
and Ha-selection functions naturally recover galaxies of 
different dust contents; if we translate these limits to 
"true" SFRs for the respective samples, we obtain limits 
of 2.6 and 6.0M Q yr -1 for the UV and Ha, respectively. 
These limits differ by a factor of over 2 in SFR, but still 
are not able to account for the differing populations of 
galaxies that survive the respective selection functions - 
were the dustier galaxies that are selected by Ha able 
to enter the UV-selected catalogues, the increased av- 
erage dust content would bring these values even closer 
together. We also argue that to some extent, the over- 
all shape of the UV and Ha LFs must be governed by 
the same physical processes and, regardless of the exact 
dust content, selecting galaxies brighter than a certain 
fraction of the characteristic luminosity should recover 
objects with similar underlying SFRs. Ultimately this 
argument is backed up in § 12. 31 when we find very similar 
UV- and Ha-derived SFRs in the local universe, and by 
th e very similar SFR densiti es derived by the tw o tracers 



tn e very similar bi 1 it aensitics derived by the two 
in IReddv etUI (|2008f ) and lHaves et all (|2010b[ ). Natu- 
rally by cutting both LFs at the same fraction of L*, we 
recover similar fractions of the luminosity density com- 
pared with integration to zero (70 %). 

For Lya, the situation is more complicated still: cut- 
ting at the same intrinsic SFR would mean that we do 
not include Lya emission at lower luminosities. This is 
now not simply a matter of dust attenuation but also 
includes radiation transport effects. Since we expect the 



line to be systematically weakened, applying a cut at the 
corresponding SFR to that of Ha or the UV would cause 
us to miss much of this light. The best way to proceed, 
therefore, is to adopt the same philosoph y as above, and 
adopt 0.04L^| 3 a . By selecting the LF of iGronwall et al.1 
l|2007h . we obtain a lower limit of 1.75 x 10 erg s . 
Should f^f c a —l. this would correspond to an SFR of just 
0.15M Q yr" 1 . However, in lHaves et al.l (|2010af) we de- 
termined a volumetric fl* a of just 5 %, and scaling this 
SFR up by a factor of 20 brings it to 2.9M Q yr -1 , almost 
perfectly into line with the UV-derived 2.6M© yr^ 1 dis- 
cussed above. Naturally, this integration from 0.04L* 
again includes ss 70% of the total luminosity density 
(compared with integrating from zero). 

In summary, selecting the optimal integration limits 
is a non-trivial process, yet we argue that by adopting 
these limits we should be selecting very similar samples of 
galaxies, at least with respect to their unobscured SFR. 
The lower limits are 4.36 x 10 27 erg s _1 Hz" 1 (UV), 4.6 x 
10 41 erg s" 1 (Ha), and 1.75 x 10 41 erg s" 1 (Lya). We 
have insured that these limits include the bulk of the 
luminosity density but are not dominated in uncertainty 
by extrapolation in the faint end, although we have also 
confirmed that integration to zero in fact has only very 
minor effects on the final measurements of / e L s J a . 

2.3. Compilation of the samples 

All of the assembled data and the derived f^. a mea- 
surements are summarized in Table [T] and Figure [T] The 
measurements of Eb-v relevant to each of the Ha or UV 
measurements are derived from data in the same publica- 
tion as the Ha or UV LF data themselves (with one ex- 
ception, which is discussed in the following paragraph). 
In this subsection we provide the necessary motivation 
for our choices and comments on the various samples. 

No instrumentation can perform a Lya-selected survey 
in the very nearby universe so we begi n at z ~ 0.2 — 0.4 
with the L ya LFs pres ented in both iDeharveng et al.l 
(2008) and ICowie et al.l (|2010t l. At these redshifts Ha 
LFs are available, and therefore we pr oceed using Equa- 
tion p We adopt the Ha LF of iTresse fc Maddoxl 
(1998), and correct it for dust attenuation by applying 
the 1 magnitude of extinct ion that is repre sentative of 
local Ha-selected galaxies ()Kennicuttlll992|) . For secu- 
rity and consistency with higher redshift measurements , 
we als o examine the z = 0.3 UV LFs of lArnouts et al.l 
(2005) , which we corr ect for dust using the method 
of iMeurer et al.1 (1 19991) a nd the /3 slope measured by 



ISchiminovich et al.l ( 20051 ) in the same sample, finding 
extremely consistent numbers. 

Beyond the very nearby universe, no further Lya infor- 
mation is available before z = 2, wh ere we adopt our ow n 
measurement of ft/ c a = 5.3 ± 3.8 % (jHaves et al.ll2010al) . 
based upon Ha and individually estimated Eb-v- 

It is already at this juncture in redshift that we lose 
the possibility to use Ha, and therefore we proceed us- 
ing published UV LFs and Equation pi Our next step is 
to take the Lya LF of ICassata et al.T(j2010l . (z) = 2.5) 
which we contrast against the dust-corrected /?l uv of 
IReddv et all (|2008l (z) = 2.3). For thi s, and all sub- 
sequent points from ICassata et all ((2010h . we adopt the 
values of L± that are uncorrected for IGM attenuation. 
It is reassuring that the measurements at z = 2.2 and 
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z w 2.5 (which are based upon Ha and UV, respectively) 
give very consiste nt numbers. Furth ermore, in a very re- 
cent submission (| Blanc et al.l I2010D an additional Lya 
LF has been presented at 1.9 < z < 2.8, the integrated 
Lya luminosity density from which differs from our own 
result by « 25 %. 

We then continue with the iReddv et all (|2008l) UV 
data at (z) — 3.05, which we use to compu t e fig " for 
the z = 3.1 Lya s amples of iGronwall et al.l (|2007[ ) and 
IQuchi et all d200l . 

At z ~ 4 we have av ailable Lya LF s from IQuchi et al.l 
(|2008l z = 3.7) a nd ICassata et all (|2010T , z = 3.9), 
and UV LFs from iBouwens et al.l (|2007L (z) = 3.8). 
We also use the z = 4.5 and 4.86 Lya LF p oints from 
iDawson et all (|2007j ) and lShiova et al.1 (|2009l ). which we 
nor malize by the dust- corrected UV point at z = 4.7 
from IQuchi et~aLl (l200l . 

The next redshift to examine is the popular z w 5.7 
Lya window. Here w e adopt the UV datap oint from the 
i— dropout sam ple of IBouwens et al.l ()2007l . (z) = 5.9), 
and the Lya LF lOuchi et all (120081. (z) = 5.7), which i s 
in good agr e ement with those ofjS himasa ku et al.1 {2006), 
lAiiki et al.1 (|2006l ). and iTapken et al l (120061). W e also 
add the highest redshift LF from ICassata etaTl pOloft 
at (z) = 5.65. 

Finally we asse mble a few z > 6 sam ples. We adopt the 
z = 6.5 p oint fromlOuchi et al.l (120101 . which includes the 
sa mple of [K ashika wa et all 12006). and the measurement 
of live et alj (l2006f) at z = 7.0, which has also been com- 
lled in lOta et al.l (|2008[ ). Here we adopt IBouwens et al.l 
2010) UV measurement at (z) = 6.8 for comparison. It 
should be noted that at this reds hift the dust-co r rected 
and uncorrected measurements of IBouwens et al.1 (|2010f ) 
converg e. We adopt the mo st optimistic estimate at z = 
7.7 from lHibon et al.1 (120101) . for which we interpolate be- 
tween th e IBouwens et all (l20Toh z = 6.8 and 8.2 UV dat- 
apoints. iHibon et al.l (|2010D present Schechter parame- 
ters for four Lya LFs, based upon various assumptions 
about the rate of contamination by lower redshift galax- 
ies. By assuming all of their candidates are real (their 
sample a) we find a Lya escape fraction of (33.5I33J) %. 
We also briefly examine their subsample b, in which only 
four of the seven objects are real. For all of their sub- 
samples, the numbers are insufficient to provide meaning- 
ful errors on the luminosity density and by our standard 
error procedure we derive / e L s ^ a =(22.2l^° 2 7 ) % Further, 
it should be noted that in the IHibon et al.1 (|2010f ) sample, 
the lower limits obtained on the Lya equivalent width 
are in the range 6-15A, with their continuum-detected 
object showing Wh ya =13k. Thus, at an acceptable con- 
fidence limit, none of their seven objects would actually 
survive the canonical Wh ya cut of 20A that is typically 
employed in narrowband surveys. Including these data is 
therefore not straightforward, but in order to treat them 
as consistently as possible with the lower redshift points, 
we have to set the z w 7.7 Lya escape fraction to zero, 
but adopted a characteristic error of 50.6 % as derived 
from their most optimistic sample. We note that this 
limit is likely extremely high. 

All of our measurements of f^ a are listed in Table Q] 
and shown graphically in Figure [1] which is the main 
result of this paper. 



2.4. Consistency (and inconsistency) between groups 

It should always be borne in mind that we are compil- 
ing results from different survey teams, who may adopt 
different techniques for data reduction and photome- 
try, derivation of the luminosity fu nctions, and incom- 
pleten ess corrections. For example, iMalhotra fc Rhoadsl 
(|2004h find reasonable agreement at z ~ 5.7 between 
the na rr owband-sele c ted L ya LFs of lRhoads fc Malhotral 
(|2001h : lAiiki et al.1 (120041) and the lensing-based sur- 
vey o f ISantos et al.l (120041). However, the z = 5.7 
LF of IShimasaku et al.l (|2006f ) , on w hich the study of 
iKashikawa et al.l ( 20061) is based (see § I5.2|) , find a strong 
disagreement at th e faint end between their ow n LF and 
the compilation of Malhotra & Rhoads (2004) . As com- 
mented by IShimasaku et al.l (|2006l ) the likely cause for 
this discrepancy lies in (a) the lack of incompleteness 
corrections, which are unmentioned in any of the 2004 
articles, (b) the differences in equivalent-width based se- 
lection crite r ia, an d (c) the large cosmic variance which 
IQuchi et al.l (|2008l ) noted can be of factors of w 2 in fields 
as large as 1 square degree. Our results are sensitive to 
all of these considerations. 

It is only now that sufficiently large samples of Lya- 
emitting galaxies are presented in the literature for this 
study to be undertaken, and we are now fortunate that 
a good fraction of our data must contain internal self- 
consistency. For example, three of our data-points at 
at z = 3.1, 3.7, 5.7) are drawn from a single paper 
Ouchi et a l. 2008) in which the methodologies must be 
internally consistent, and the basic trend can be seen 
in th ese data alone. A fourth point at z — 6.6 comes 
from IQuchi et al.l (|2010T l where similar self-consistency 
is to be expected. I n the same fashion, three further 
points are taken from ICassata et all (j2010| ) where inter- 
nally the same methodology must have been adopted at 
each redshift. It is certainly encouraging that, for ex- 
ample at z = 5-7 the measurements of /iff" b ased upon 
ICassata et all (|2010f ) and IQuchi et all (|2008f) are prac- 
tically indistinguishable, despite the fact that they are 
based upon completely different methods: blind spec- 
troscopy and narrowband i maging, respec ti vely. The 
z = 2.2 and 2.5 p oints of ICassata et al.l (|2010f ) and 
I Haves et al.l (|2010a|) are similarly indistinguishable, as 
(and also robust against the same fundamental method- 
ological difference of blind spectrosc opy vs narrowband 
ima ging), are the z = 3.1 points of IQuchi et al.l ()2008l ) 
and IGronwall et all (|2007l ) (both narrowband imaging). 

Any study of the galaxy population benefits by target- 
ing spatially disconnected, independent pointings in or- 
der to beat down cosmic variance. By adopting the stud- 
ies of various authors pointed all over the extra-galactic 
sky, this study is able to benefit from the inclusion of a 
large number of independent fields. 



3. GENERAL RESULTS 

3.1. The evolution of flg a 

Figure[T]reveals a general and significant trend for f^/ c a 
to increase with increasing redshift. Beginning in the 
very local universe we see /^°~ 0.01 or lower for nearby 
star-forming objects. This increases to around « 5 — 10 % 
by redshift of w 3-4, and further to « 30 - 40 % 
by redshift 6. In order to quantify this trend we fit 
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TABLE 1 

Lyman-alpha escape fractions with redshift 





Lyct quantities 




1 


Intrinsic quantities 




| Derived results 


z 


Ref 


P* 


z 


Ref E B -v 


P* 


/S° [ % ] Comment 


(1) 


(2) 


(3) 


(4) 


(5) (6) 


(7) 


(8) (9) 



Estimates based upon Lya and Ha luminosity functions 



0.2-0.35 De 08 
0.2-0.4 Co 10 
2.2 Ha 10 



(3.79 ± 1.69) x 10~ 4 
(8.33 ± 2.60) x 10~ 5 



0.2-0.35 
0.2-0.35 
2.2 



TM 98 
TM 98 
Ha 10 



0.33 
0.33 
0.22 



(0.0303 ± 0.017) 
(0.0303 ± 0.017) 



(1.25 ±0.90) 
(0.275 ±0.18) 
(5.3 ±3.8) 



1 mag at Ha 
1 mag at Ha 
Multi dimensional M.C. 



Estimates based upon Lya and UV luminosity functions 



2.5 


Ca 10 


(7.08 ±0.81) x 10~ 


-3 


(2.3) 


Re 08 


0.15 


(0.201 ± 0.022) 


(3.51 ±0.56) 




3.1 


Gr 07 


(8.50 ± 5.32) x 10" 


-3 


(3.05) 


Re 08 


0.14 


(0.116 ± 0.017) 


(7.33 ±4.71) 




3.1 


Ou 08 


(5.54 ±2.91) x 10" 


-3 


(3.05) 


Re 08 


0.14 


(0.116 ± 0.017) 


(4.78 ± 2.61) 




3.7 


Ou 08 


(4.78 ± 1.14) x 10" 


-3 


(3.8) 


Bo 09 


0.14 


(0.089 ± 0.011) 


(5.36 ± 1.43) 




3.8 


Ca 10 


(8.71 ± 1.00) x 10- 


-3 


(3.8) 


Bo 09 


0.14 


(0.089 ± 0.011) 


(9.77 ± 1.64) 




4.5 


Da 07 


(3.22 ± 1.25) x 10- 


-3 


(4.7) 


Ou 04 


0.075 


(0.025 ± 0.011) 


(12.6 ± 7.17) 




4.86 


Sh 09 


(2.35 ±3.17) x 10" 


-3 


(4.7) 


Ou 04 


0.075 


(0.025 ± 0.011) 


(9.24 ± 13.0) 




5.65 


Ca 10 


(8.53 ±3.44) x 10- 


-3 


(5.9) 


Bo 09 


0.029 


(0.022 ± 0.005) 


(38.1 ± 17.2) 




5.7 


Ou 08 


(6.76 ±4.77) x 10" 


-3 


(5.9) 


Bo 09 


0.029 


(0.022 ± 0.005) 


(30.2 ± 22.2) 




6.6 


Ou 10 


(4.73 ± 1.24) x 10- 


-3 


6.5 


Bo 07 


0.012 


(0.016 ± 0.008) 


(30.0 ± 17.8) 


UV Interpolated 


7.0 


iy 06 


(1.07 ± 1.16) x 10" 
(0+g 8 ' 5 ) x 10- 


-3 


7.0 


Bo 09 


0.010 


(0.012 ± 0.008) 


(8.96 ± 11.5) 


UV Interpolated 


7.7 


Hi 10 


-3 


7.7 


Bo 10 


0.0 


(0.005 ± 0.002) 


(0 +so. 6) 


UV Interpolated 


Note. 


— For the 


Ha-bascd estimates, i 


ve use the 


integrated 


luminosity densities directly; SFRD measurements are presented 



for homogeneity wit h the UV estimates. p+ units of are yr 1 M pc d and Eb-v is in magn itudes. The references are ex - 

panded as: Bo 09=IBquwens et al.1 d2009aT): Ca 10=|Cassata et all 11201011: Co HMCowie et al.l <201Cfl: Da CTMDawson et all ||2007I 



ll2010aft : Hi 09=IHibon et al.1 f 20 1 Of) : Iv 08=IIve et at] 



De 08 = lDeharveng etal . (2008); Gr 07=Gronw all et al 
( 120061) ; Ou 04= IOuchi et al.1 112004ft ; Ou 8= IOuchi et al.1 _ 
( 120081) ; TM 98 HTresse &: Maddoxl i 1998ft . References for E B -v measurements are the same as for the intrin sic star-formation rate density 
(I.e. that listed in the 5th column) with the exception of the (z) — 0.3 points in which Eb-v is adopted fromjKcnnicutt ( 1992ft . 



(2007); Ha 10= Haves et al. 

2008ft ; Ou lO MOuchi et al.l i hoiOftT Sh 09 MShiova et al.l ( 120091) ; Re 08 MReddv et al.l 



an analytical function to these data-points, choosing a 
power-law of the form fcZ. a ( z ) = C • (1 + z)^ - we obtain 
coefficients of C = (1.67±g;||) x lO^ 3 ;^ = (2.57±g;il). 
Note that we do not include any z > 6 points in our 
fit since previous studies suggest that it is around this 
redshift that an appreciable fraction of the intergalac- 
tic hydrogen becomes neutral, and may in principle af- 
fect the Lya LF. For more discussion on this see § 15.21 
To insure that the fit is not biased by the presence of 
two z « 0.3 points that lie around 8 Gyr from z w 2, 
we repeat the fit after excluding these points, finding 
C = (4.79± 'J|) x 10- 4 ;£ = (3.38±g;$). Clearly the fit 
is affected by these points, but their exclusion actually 
results in a more rapid evolution with redshift. 

Beyond redshift 6 the apparent trend begins to break 
but it is initially very slow. Over the redshift interval 
of 5.7 to 6.5, fe/ c a stabilizes, but decreases aga in to just 
ps 10 % at z = 7. The redshift 7 point from live et al.l 
(2006) is confirmed, whe reas none of the sa mple of red- 
shift 7.7 candidates from lHibon et all (|2010T ) have confir- 
mations by spectroscopy, and this upper errorbar must 
be regarded as an optimistic upper limit. 

Finally, we perform a simple experiment with the best- 
fit relationship to the fe/ c a —z trend, and extrapolate to 
estimate the redshift at which f^. a reaches unity. This 
would carry the implication that the ISM of the aver- 
age galaxy has become effectively devoid of dust, and 
since dust is a byproduct of the star-formation process, 
must also correspond to a time of approximately primeval 
star formation. It is interesting, therefore, that we find 
/osc a = l at z = ll.ll^g, which is consistent with the 
redshift of the instantaneous reionization of the Universe 
based upon W-MAP data [z = 11 ± 1.4; iDunklev et all 



120091) . 

3.2. Comparison with the literature 

Naturally this is not the first time that f^ a has been 
estimated and several other studies based on a wide ar- 
ray of methods have attempted to pin down the same 
quantity at different redshifts. 

For example, at redshifts of 5.7 and 6.5, we compute 
/esc™ 01 around 40 % and 30 %, respectively. Based 
upon the fitting of spectral energy dist r ibutio ns (SED) 
to stacked broadband fluxes, lOno et ail (|2010D estimate 
/Ly«= (36l^) % and (4tJ 8 8 °) at the same redshifts. Al- 
though derived from an interesting approach, the uncer- 
tainties are still too large to p rovide a useful comparison. 

Like us . INagamine et ah (|2008|) compared observed 
Lya LFs (|Ouchi et al.l 12006 , in this case, which we also 
use) with intrinsic estimates, having derived this intrin- 
sic LF from smoothed particle hydrodynamical (SPH) 
models of galaxy formation. They adopt two methods 
of scaling the intrinsic to the observed LFs, the first of 
which they call 'escape fraction', which is a scaling to 
the datapoints along the luminosity axis, and assumes 
all galaxies have the same f^ a . This method finds 
/esc a= 10 % at 2 = 3, which is certainly consistent with 
our es timates based on the z — 3.1 LF of iGronwall et all 
(|2007l ) and sim ilar to, but s l ightly higher than our esti- 
mate based on lOuchi et aLI (|2008l ). At z = 6 however, 
INagamine et al.l (|2008[ f require an escape fraction of just 
15 % which is lower than our estimates of 30 — 40 %, and 
discrepant with our e stimates at around the 2<r level. 
INagamine et al.l (|2008l ) also test a 'duty-cycle' scenario 
(an LF scaling along the $ axis) in which only a frac- 
tion of the SPH galaxies are 'on' as Lya-emitters but 
emit 100 % of their Lya photons. Note that in these two 
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Fig. 1. — The redshift evolution of /, 
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Publication cod es are listed in the footnote to Table [T] 



3.1 and 5.7 points have been 



artificially shifted by Az = 0.08 for clarity. The point from jHibon ct al. 2010]) takes, according to our definition, a value of zero. It is 
therefore displayed at a value of 0.002 to permit visualization on a logged axis. The solid red line shows the best fitting power-law to points 
between redshift and 6, which takes an index of £ = 2.6 and is clearly a good representation of the observed points over this redshift 
range. It intersects with the /<J J s ? a =l line (dotted) at redshift 11.1. 

extreme scenarios, there is no requirem ent for the inte- 
gral oy er the scaled LF to be equivalent. FNagamin e et al.l 
(2008) present duty cycles of 0.07 and 0.2 at z = 3 and 6, 
respectively. However, before they compute these scal- 
ings the observed LFs are shifted along the luminosity 
axis by IGM attenuation factors of 0.82 (z = 3) and 0.52 
(z = 6), which also need to be applied for a compari- 
son with our estimate. Thus in the duty cycle scenario, 
the volumet ric escape fractions that one would infer from 
the study of lNagamine et all ()2008l ) are 6 % at z — 3 and 
10 % at z = 6. Again this agrees very well with our mea- 
surement at z » 3 but compared with our estimates at 
z — 6 is an underestimate of around the same magnitude 
as their escape fraction method. 

In contrast, using similar SPH galaxy formation mod- 
els but modified prescriptions for Lya production and 
transmission, as w ell as a different reionization history, 
iDaval et al.l ((2009) find Lya escape fractions of 30 % at 
both z = 5.7 and 6.5, which corresponds exactly with our 
measurements. Similar values of /<^ Q ~ 2 3-33% have 
also b een obtained in the follow-up work of IDaval et al.l 
(2010), although they include also an IGM transmission 
of T a — 0.48. Throughout this paper we have made sure 
not to apply any IGM correction, since the value of T a 
remains poorly constrained, even theoretically, and from 
an observational perspective there is no strong evidence 
for exactly how c lose the IGM comes to a narrow Lya 
line. As with the iMadaul (|1995[ ) prescription, it is likely 
that this IGM transmission is too low when considering 
lines that are systematically redshifted by the kinemat- 
ics of the ISM, which would drive up these theoretical 
estimates of the Lya escape fraction. 

Adopting a similar m ethod of LF scaling by luminosity, 
iLe Delliou et al.1 ([20051 ) found that an escape fraction of 
2 % was sufficient to match observed Lya LFs with their 
predictions based upon semi-analytical models between 
z = 2 and 6, with the same machinery able to predict the 



clustering properties of Lya emitters ([Orsi et al.l [2008). 
This is at the lower end of being consistent with our z = 3 
measurements, and should the same escape fraction hold 
at z = 0.3, would also be consistent wi th our estimates 
in the nearby universe. However, the ILe Delliou et all 
( 2005) escape fraction is highly inconsistent with our es- 
timates at higher redshi ft. These sem i -analy tical models, 
using the prescription of lBaugh et al.l (|2005l ) , categorized 
star-formation as occurring in two discrete modes, with 
a normal Salpeter IMF (a = —1.35) assigned to quies- 
cent star-formation and a flat IMF (a = 0) for bursting 
systems. This flat IMF increases the ionizing photon 
production at a given SFR by a factor of ten and was 
implemented as a requirement in order to reproduce the 
population of su b-mm selected ga l axies at z > 2. How- 
ever as noted bv ILe Delliou et all (|2006l ). the fraction of 
total star-formation that occurs in bursts increases from 
5 % at z = to over 80 % at z — 6, and thus their 
model implies that by the z ~ 5.7 points, effectively 
all stars are formed in environments where ionizing pho- 
tons are greatly over-produced compared to the present 
day. However, should this requirement of the flat IMF be 
removed and Salpeter applied throughout, the intrinsic 
rate of production of ionizing photons would be decreased 
by a factor of 3 at z = 3.1 where the star- formation is 
shared evenly between bursting and quiescent systems. 
This would bring the estimate to 11 % at this red- 
shift. At z = 6, /^"=16 % would be found by replacing 
the flat IMF with Salpeter. Th ese numbers are i ndeed 
very similar to the SPH models of Na gamine et al.1 (2008) 
but inconsistent with those of lDaval et al.l ( 20091 ) and our 
own estimates based upon observation. It is interesting 
to point out, however, that the IMF assumption has lit- 
tle effect on the z ~ 0.3 points where, in their model, the 
quiescent mode of star-formation dominates. 
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3.3. Possible physical explanations 

The evolution in measured f^. a is substantial, cov- 
ering approximately two orders of magnitude, and no 
doubt holds vital information about the physical nature 
of galaxies at various cosmic epochs. As we will show in 
§ [5j the most likely explanation for this evolution is the 
decrease of the average dust content of galaxies. However 
from a physical perspective many effects may enter. For 
example, galaxies may also contain less neutral hydro- 
gen to scatter photons, show faster outflows, or become 
more clumpy. The inferred increase may alternatively 
be mimicked by galaxies becoming younger on average, 
having low and decreasing mctallicities, or forming stars 
with IMFs that become more biased in favor of mas- 
sive, ionizing stars. On the other hand, the scattering of 
Lya photons by a neutral IGM, and the general leakage 
of ionizing photons (LyC) are expected to increase with 
increasing redshift, and would both serve to lower the 
perceived Lya escape fraction (although the "true" f^f c a 
of galaxies, i.e. before the IGM, would not be affected). 

Regrettably we are not able to measure any of these 
quantities directly from this compilation of data. We 
have, however, assembled data that show a number of 
trends with redshift: the Lya and UV luminosity densi- 
ties and the dust contents. These we have combined to 
show how evolves, yet in order to extract the maxi- 
mum of information from these, we need to examine an- 
other possible trend: how correlates with dust con- 
tent. Thus we delay a detailed discussion of what drives 
the fl& a -z trend until § and now proceed to discuss 
the effects of radiation transport and dust absorption. 
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Fig. 2. — Literature co mpilation of / e L S c° v s Eg_y. The codings 
in the legend arc: Ha 10= Haves et al. (2010a); Ko 10=Kornci ct al. 
(2010). Solid circles from Haves ct al. (2010a) are six objects for 
which we have detections in both Lya and Ha. Caret down markers 
are Ha emitters that were undetected in Lya and hence presented 
as upper limits, while caret up markers are Lya galaxies for which 
Ha lies below the detection limit and are hence presented as lower 
limits. Errorbars are removed from the plot to aid readability, 
but the average errors from the common detections of lHaves et al.l 
(2010a) are shown by the singular black point with er rorbars. For 
further information the reader is referred to Figure 3 of lHaves et al.l 
(2010a). The red lines show various conversions between the ob- 

served stellar E^—v and / C L S J° . The dotted line shows the standard 
ICalzett i ct al. (2000) prescriptio n, the dashed line s hows the 1 di- 
mensional fit to the data from [Hayes ct al. (2010a) and the solid 
line a 2 dimensional fit described in the text. 



4. THE Lya ESCAPE FRACTION AND ITS DEPENDENCIES 

That Lya photons undergo a complex radiation trans- 
port, in which a large number of parameters enter, is 
well-known but poorly understood from an empirical an- 
gle. Transpo rt is thoug ht to be affec ted by dust content 
(lAtek et all I2008L 120091 lHaves et al.ll2010al ). dust geom- 
etrv dScarlata et al.1 12009[). Hi content and kinematics 
(iKunth et al.lll 998: Mas -Hesse et al.l feOOS: Sh aplev et al.1 



l2003tlTapken et al]|2007 ) , and g e ometry /neutral- i onized 
gas topolog y (iNeufeldl Il99ll; iGiayalisco et ail 119961: 
Hansen fc Obi 120061: iFinkelstein et al.l 120081 120091) . Un- 
fortunately, Hi masses remain impossible to measure di- 
rectly beyond the very local universe. Kinematic mea- 
surements of the neutral ISM can be obtained at high- 
rcdshift, but require deep absorption line spectroscopy 
against the vanishing continuum of Lya-selected galax- 
ies and thus are prohibitively expensive for large samples 
of individual galaxies. We are therefore effectively lim- 
ited, when targeting statistically meaningful samples, to 
examining Lya emission against the dust content, and 
have to infer information about the remaining quantities 
by secondary analysis. 

Significant anti-correlations between f^/ c a and Eb-v 
have been presented in four recent papers, all of which 
invoke different selection funct ions and employ differ- 
ent methods of analysis. Firstly. IVerhamme et al.l (|2008l ) 
used radiation transport modelling of spectrally resolved 
Lya features in a sample of LBGs between redshift 
2.8 and 5 to estimate both dust attenuation and f}Z- a ■ 
Based upon the Balmer line ratio (Ha/H/3), Ate k et all 



(2009) computed f^/ c a and nebular reddenings based 
upon purely nebular physics in a sample of nearby 
Lya-selected galaxies. Were Ha and H/3 observations 
available in the distant universe, this method would 
be the ideal one by which to proceed. More recently, 
iKornei et al.l (|2010n performed a similar experiment in 
a sample of redshift ~ 3 Lya-emitting LBGs, in which 
dust attenuation and intrinsic Lya luminosities were es- 
timated from modelling of the SED. Finally in sample of 
redshift 2 Lya- and Ha-selected galaxies, we also used 
SED modeling to estimate Eb-v but estimated the in- 
trinsic Ly a production from the dust-corrected Ha lumi- 
nosity paves et al.ll2010af ) . 

In Figure [2] we show a comp il ation of t he and 
Eb-v points from IKornei et all (|2010T ) and lHaves et al.l 
(2010a). Here we adopt only these two data-sets since 
they involve similar computations of Eb~v but include 
Lya, Ha, and UV selection and should be broadly rep- 
resentative of the general galaxy populations under con- 
sideration in this paper. These two studies both per- 
form full SEP fits, but use them in different ways, with 
IKornei et all (|20iq) requiring the intrinsic ionizing pho- 
ton budget to estimate and lHaves et alj (|2010al) 
using only the Eb-v e stimate to correct H a for the dust 
attenuation. Thus the IKornei et al.l (|2010D points are in 
principle expected to be more sensitive to the standard 
set of assumptions in population synthesis (IMF, stellar 
atmosphere models, etc). However a substantial overlap 
between the two populations is clear in Figure^ with the 
two populations occupying a very similar region of the 
fesc a ~^B~v plane (the fact that we find more galaxies 
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at higher Eb-v is due to the fact we find redder galax- 
ies by Ha selection than is possible using the UV-biased 
Lyman-break criterion). 

The dotted line shows the dust attenuation prescrip- 
tion of Calze tti et al.l (|2000f) which should be valid in 
the case of no Lya scattering and a simple dust screen. 
This line is described by f^g a = Kr - 4 '^-^"^ ; w h e re 
^1216 = 12- Very few points lie above this line and all 
are likely placed there by statistical scatter. Indeed, this 
line sets an approximate upper limit to the datapoints, 
which extends in the direction of lower f^. a due to radia- 
tion transport effects increasing the effective dust optical 
depth seen by Lya. 

In attempts to quantify the effects of r esonance scatter- 
ing and dust absorption the s tudies of iVerhamme et al.1 
(|2^ . lAtek~eTanj|2^09f ). and lHaves et al.l (|2010aD all fit 
linear relationships to the datapoints on the log(/^ Q )- 
Eb-v plane, assuming no a priori information about 
the dust. These studies all used a functional form of 
/Ly Q = iQ-0-4-E B -vk Lyai ^ wn ere k Lya (the single free pa- 
rameter of the fit) is an effective extinction co-efficient for 
Lya, and thus includes both sc attering and abso r ption . 
Bot h at high-z, the stu dies of rVcrham me et al.l (2008) 
and lHaves et al.l (|2010al) found effectively the same value 
of ki.vn— 17.8, which runs significantly steeper than the 
ICalzetti et al.l (|2000f) relationship as Lya photons are 
preferentially attenuated. This is shown by the dashed 
line in Figure [5J 

These formalisms force the fits to conform to fe^ a =l 
at Eb-v=®, and technically it is true that if there is ex- 
actly zero dust, Lya photons cannot be absorbed by dust. 
However, the very presence of Lya photons implies that 
star-formation must be occurring and, after just ~ 3 Myr 
of star-formation, dust produced in supernovae would be 
returned to the ISM and the optical color excess ceases to 
be a good proxy for dust. It is well-known that Lya can 
be strongly suppress ed even when miniscule amounts of 
dust are present (e.g. lHartmann et al.lll984HKunth et al.l 

[1991 iThuan fc Izotovl Il997l: IQstlin et al.l 12001 ) and as 

Figure [2] shows some galaxies have fcZ a= ^®% with no 
measurable UV attenuation. Indeed, many star-forming 
galaxies show little or no attenuation in front of their 
ionizing clusters but substantially attenuated nebular 
regions. This is the origin of the factor of 2.2 diffcr- 
ence between stellar and nebular measurements of Eb-v 
(|Calzetti et al.l [2000h . but at a very low UV stellar at- 
tenuation of Eb-v~ applying a factor of two is not 
meaningful and nebular lines in general - and Lya in 
particular - may be heavily attenuated. It is unfortu- 
nate that at high-z the UV continuum is our only proxy 
for the dust content as we indeed expect to be surveying 
redshifts at which the stellar a ttenuation indeed falls to 
~ fe.g. lBouwens et al.|[2009al) . 

To account for these factors we now proceed to 
relax the requirement of the fit passing through 
(Eb-v ,/jff a )=(0,l) and re-fit the combined datasets of 
IKornei et all pOlof ) and lHaves et al.l (|2010al) using the 
following expression 
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(4) 



This expression takes the same form as the standard 
dust-screen prescription, with coefficient fcLya, but adds 
the additional free parameter of CL ya , the factor by 



which f^/ c a is scaled down. As in lHaves et all (|2010aJ) 
we use Schmidt's binned linear regression algorithm 
(jlsobe et al.l ll986). since it permits the combination of 
data-points and limits in both directions. For fcL ya we 
obtain a value of 13.8, which is much more sim i lar to 
the value of 12.0 obtained from ICalzetti et"aT1 (|2000l ) 
at the wavelength of Lya. However, we also obtain 
CLya=0.445, indicating we expect f^/ c a to be around 
50 %, even when there is no measurable dust attenuation 
on the stellar continuum. This is in fact a more plausible 
scenario since the effect of scattering by neutral hydrogen 
is not expected to depend on the dust content itself. This 
fit is show n by the soli d red l ine in Figure [21 A g ain the 
points of IKornei et all (|2010l) and lHaves et all (|2010aft 
are subject to different assumptions that enter the pop- 
ulation synthesis. However, for the reasons outlined pre- 
viously in this subsection and the similarity between the 
distributions, we do not expect these quantities to be 
strongly subject to these assumptions. 

It is not necessarily straightforward to define a 
goodness-of-fit measurement to compare the quality of 
the three fits, given the large number of upper- and lower- 
limits in this dataset. Thus we define our own normalized 
r.m.s. statistic (rms n ), as: 
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where /j meas is the i measured Lya escape fraction, 
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EBV 



IS 



the i th 



Lya escape fraction predicted from 



Eb-v, and N is the number of data-points. However 
in order to treat the limits, we permit a point to con- 
tribute to the summation only if that limit is violated. 
We appreciate that this is a non-standard statistic, but 
it does enable a quantitative measure of the goodness-of- 
fit that is philosophically not too far removed from more 
commonplace stat istics. Adopting the f l^-En-v rela- 
tions derive d from ICalzetti et aLl ( 2000f ) . the one param- 
eter fit from lHaves et al.l ()2010aD and the two parameter 
fit from this work, we compute rms n — 1.85, 1.02, and 
0.66, respectively. 

We have now assembled information about three 
trends: the observed redshift evolution of f^. a ; the ob- 
served redshift evolution of the dust content of galaxies; 
and the observed relationship between f^/ c a and dust 
content. We will next show that we are able to syn- 
thesize these points to infer some general trends in the 
evolution of galaxies. 
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5. ON THE EVOLUTION OF 

5.1. Redshifts 0-6: the upwardly evolving escape 
fraction and the properties of galaxies 

5.1.1. The evolving dust content of galaxies 

We showed in the previous section that of in- 

dividual galaxies is anti-correlated with the measured 
Eb-v (Figure [2]). Given that the typical Eb-v evolves 
with redshift (see Table [T]), we may indeed expect a pos- 
itive correlation between and redshift. This is ex- 
actly what Figure Q] shows, where it is clear that the 
Lya escape fraction increases smoothly and monotoni- 
cally out to z ~ 6. Thus it appears that this increase in 
/esc" i s the result of the dust content of the star-forming 
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redshift 

Fig. 3. — Same as Figurefl] but zoomed onto the relevant region. 
The red lines show the Lya escape fractions that would be pre- 
dicted based upon the values of Eb—v that have been measured in 
the respective Ha and UV samples (listed in TableflJ, and using the 

various conversions between measured Eb—v an d fcsc* described 
in the text. The d otted line represents the dust attenuation law of 
ICalzett i ct al. (2000), the dashed line the 1 dimensional empirical 
fit to the data of Haves ct al. (2010a), and the solid line a 2 di- 
mensional fit to the data described in § [J] Using the 2 dimensional 
fit, a remarkably good agreement is seen between observations and 
prediction between redshifts and 6.5. 

galaxy population decreasing with redshift. We now take 
the measured values of Eb-v from the various samples 
(listed in Table [1} , and use them to compute the f^. a 
that would be expected, from the three conversions be- 
tween Eb-v and f}£ a discussed in the previous section 
[Calzetti et al.l (120001) an em pirical fit with one free pa- 
rameter ( Haves et al.ll2010al ). and an empirical fit with 
two free parameters] . We show the measured escape frac- 
tions together with these predictions in Figure O 

We first d is cuss the predictions based upon the 
ICalzetti et al.l (|2000L red dotted line), which is clearly 
discrepant with the observations at around the 3cr level 
at every redshift. Obviously this is to be expected since 
Lya photons resonantly scatter and it is unlikely that 
the dust is distrib uted in a uniform sc reen. The one di- 
mensional fit from lHaves et al.l (|2010a[ ) offers substantial 
improvement and is able to describe the observations be- 
tween redshifts and 4. This reasoning is circular for 
the redshift 2 points where the f^^-Es-v relationship 
was derived, but we stress the tautology is present only 
at this redshift. This relationship is not able to explain 
any of the datapoints at redshift above 4, where it sys- 
tematically over-predicts the Lya escape fraction. 

As redshift increases the dust content of galaxies is 
clearly shown to change and, could we plot Figure [5] at 
redshifts higher than 3, we could expect galaxies to clus- 
ter successively fu rther toward s the up per left corner of 
the plot. Since the lHaves etall (|2010a| ) f^ c a -E B -v fit is 
forced through the {Eb-v Jci c a )— (0,1) coordinate and a 
high value of k^ ya is found, the predicted escape fraction 
evolves very quickly with redshift. Indeed, these predic- 
tions evolve much faster than the data, as is forced 
for unphysical reasons towards unity. 

When we introduce the new f^/ c a -EB-v fit with two 
free parameters and allow C\, ya ^ 1, the agreement be- 
tween the measured and observed Lya escape fractions 
is striking: it agrees with essentially every datapoint, 



within the errorbars, between redshift and 6.6. We 
should point out that it is not clear that the use of the 
average Eb-v for a sample should by necessity repro- 
duce the volumetric escape fraction. Due to variations of 
the dust contents and ISM of individual galaxies, and the 
associated impact upon the transfer of Lya and the selec- 
tion of galaxies, it is possible that the average f^/ c a could 
have been skewed substantially from the data-points. In- 
deed, close examination of the fl"Z a -En y rela tionship 
(CLya=l; fcL yQ =17.8) from lHaves et al.l (|2010af ) reveals 
that it does not perfectly intersect the center of the f^/ c a 
datapoint (z — 2.2 point in Figure |3]) from the same sur- 
vey, despite /<J^ a , average Eb-v, and the coefficients of 
the f^/ c a -Es-v relationship all having been derived en- 
tirely from this one dataset. This most likely results from 
the weighting across the population from which the aver- 
age Eb-v is computed (the representative E B -v is not 
an average weighted by the intrinsic Lya luminosity) , ex- 
actly the effect under discussion. However, the fact that 
such tight agreement is seen between the observational 
estimates and those derived from our fit suggests that 
such a bias in the selection of the populations is not at 
play here. 

Again we stress that the relationship we derived be- 
tween fc/ c a and E B -v in §IU includes the effects of res- 
onance scattering, and thus in some manner the neu- 
tral gas content, its kinematics and relative geometry all 
enter the relationship, which holds even when the mea- 
sured optical color excess on the stellar continuum is zero. 
There is no reason to assume that these quantities are 
constant with redshift and we could, for example, envis- 
age situations where the gas content, feedback properties, 
or dumpiness evolve and thereby change fcLya or C\^ ya . 
However the tight agreement between our observed 
values and those computed from the f^^-Es-v rela- 
tionship provides no evidence for the evolution of these 
properties (at least if the gas content does change it does 
not take part in the Lya scattering process). The evolu- 
tion of fl* a across almost the entire observable universe 
can be explained cleanly within the confines of this sim- 
ple model, as mainly due to a dust content that evolves 
with redshift. 

5.1.2. Other effects 

We need to interpret an increase in the global 
of galaxies by a factor of ~ 4 between z = 2 and 6, and 
naturally if something were to alter the intrinsic Lya/UV 
ratio of galaxies by this factor, the evolution in f^. a 
could be mimicked. 

For example, there is evidence that the W\^ ya distribu- 
tion of galaxies changes with increasing redshift: high- 
Wlvq objects becom e rel atively more abundant (e.g. 
Gronwall et~aLl 120071 c.f. iShimasaku et al.l 120061: also 
Ouchi et al.l 120081 ) . and thus pure selection may explain 
the trend. However, the Wh ya distributions at z — 2 and 
3 suggest a maximum of ~ 20 % of the total luminosity 
density will be lost by non-selection of < W\, ya < 20 A 
galaxies, and such a selection bias can certainly not ex- 
plain the magnitude of the trend observed here. 

It may also be argued that lower metallicities or a flat- 
tening of the IMF may explain the trend. However, 
between solar and 1/50 solar metallicity the increase 
of Wh ya for constant SFR, a measure of the relative 
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Lya/UV output, is less than 50 % (jRaiter et all 120101 ). 
insufficient to explain the observed increase of f^. a . To 
explain an increase by a factor ~ 4 would require a de- 
crease of the ave rage metallicity fro m solar down to less 
than 10 -3 solar (jRaiter et al.ll2010f ). which seems highly 
unlikely. 

One would also assume that a relatively higher frac- 
tion of genuine primeval galaxies would be discovered as 
rcdshift increases, and a substantial (~ 3-fold) enhance- 
ment of Lya/UV may arise from preferential selection 
of extremely young systems (e.g. iCharlot fc Falll 119931 : 
Schacrcr 2003). To get this kind of enhancement a galaxy 
must either be observed at an age below ~ 10 Myr or, 
should an episode of star-formation occur superimposed 
atop an aged stellar population, sufficient time must have 
elapsed for that population to fade in the UV. For this 
UV fading to occur, punctuated bursts of star-formation 
would need to be separated by around the UV equilib- 
rium timescale of ~ lOOMyr. At z — 6 the Universe 
has an age of 1 Gyr and even if all star-formation were 
to occur in individual bursts, the chance of catching an 
individual galaxy at this time would be around 10 %. 
Thus, integrated over the entire galaxy population the 
application of such a sampling bias also seems quite im- 
plausible. 

We may expect at some point over this cosmic evolu- 
tion, that galaxies start to leak a substantial fraction of 
their ionizing photons {fc/ c C )- Indeed as we approach 
the middle of the epoch of reionization, the reionization 
processes itself dictates that this must be true, and we 
may expect at lower redshifts (e.g. 4-6) that a substan- 
tial population of galaxies may remain with an ISM that 
permits high /^ c . In addition, across approximately the 
same redshift domain we may expect the thickening neu- 
tral phase of the IGM to start to suppress Lya. Both of 
these effects would act to lower the perceived Lya escape 
fraction by either draining ionizing photons or scattering 
Lya. Although we are not able to tell whether these 
effects become significant at z ~ 4 — 6, if they do be- 
come important then the intrinsic Lya escape fractions 
of these galaxies will be still higher than we measur^. 

It may be argued that the measured Lya fluxes (and 
hence the Lya luminosity density) could be underesti- 
mated due to the spatial extension of Lya, and that 
some of the observed redshift trend could be d ue to 
this (e.g. ILoeb fc Rvbickil U99<l IZheng et alJI2010D . Al- 
though a somewhat larger spatial extension of Lya com- 
pared to th e UV continuum has been noted in som e 
surveys (e.g. lNilsson et al.ll2009t iFinkelstein et al.ll2010D . 
stacking analysis in other Hubble Space Telescope im- 
ages reveals the Lya emission to be spatially compact, 
with only a small fraction of the integrated luminos- 
ity lost to aperture effects ([Bond et al l 120101 ) . There- 
fore it seems very unlikely that this could lead to a 
significant underestimate of the Lya flux, which would 
mimic the apparent trend of increasing Lya escape frac- 
tion with redshift. The main reasons are the follow- 
ing. First, the photometric apertures typically used 
for the narrowband images taken from the ground are 

8 For example, assuming that half of the Lya flux is lost due to 
scattering in the IGM the "intrinsic" value of out of galaxies 

would be higher by a fa ctor 1.22 (1.92 ) at z ~ 3 (6), assuming the 
average IGM opacity of Madau (1995). 



several times larger than the FWHM of the Lya emis- 
sion and several studies apply the same me thod at sev- 
eral redshifts (e.g. between z ~ 3 and 6, lOuchi et al.l 
120081) . Second, several independent measurements using 
both imaging and spe ctroscopy reveal the same trend 
betw e en z ~ 2 and 6 (|Ouchi et all 120081 : ICassata et al.l 
2010t iStark et al.ll2010al). and also over a smaller red- 
shift range (|Reddv et aLll2008t ). Third, it is well-known 
that in individual Lya-selected systems at redshifts 2-3, 
the SFR inferred by comparing Lya and UV radiation is 
frequently found to be comparable (Guaita et alj|20lot 
INilsson et ai]|20t)9t lOuchi et all 120081) . Finally, some of 
the brightest Lya-emitting objects on the sky - where the 
order-of-magnitude fainter low surface brightness scat- 
tered emission should b ecome apparent - a lso seem to be 
spatially compact (e.g. lWestra et alll2006l ). These obser- 
vational lines of evidence all argue against an important 
loss of Lya photons related to its spatial extension. 

At z ~ 0.2 — 0.3 the Lya emitting samples have been 
carefully constructed from surveys using GALEX slit- 
less spectroscopy of NUV continuum se lected objects 
(jDeharveng et al.ll2008HCowie et alll2010D . Given its rel- 
atively low spatial resolution (~ 5") the Lya flux mea- 
surement of individual sources should not be affected by 
possible differences in the spatial extension. Further- 
more, blending affe c ts on ly 10% of the sources, accord- 
ing to iCowie et al.l (|2010l ) . Finally, comparing number 
counts of GALEX sources with/without Lya emission 
these authors have also shown that the Lya emitters 
represent only ~ 5 % of the NUV-selected continuum 
sources, a fraction significant ly lower than t he 20 -25% 
derived for z ~ 3 LBGs by IShaplev et all pOOll ). In 
other words, a low escape fraction at low-z is not only 
obtained from the ratio of the UV and Lya luminosity 
density, but also from direct inspection of NUV contin- 
uum selected objects. 

Finally, as discussed in § 12.21 our assumed limits of in- 
tegration may introduce an overall bias into the data. For 
both the Lya- and UV-selected populations, the char- 
acteristic luminosity of the LF (L+) is known to evolve 
with redshift. Thus selecting a constant lower limit at all 
redshifts may result in an artificial evolution. Firstly it 
should again be noted that our fixed lower limits apply 
to both the numerator and denominator (Lya and UV 
LFs; in Equation [1]) and to first order will cancel. Sec- 
ondly, the evolution of both Lya and UV LFs follows a 
similar pattern, starting low in the nearby universe and 
increasing rapidly to z = 2 or 3, from where they be- 
gin to decline in the direction of the highest redshifts 
(with the Lya LF declining slower than that of the UV 
in this range). Thus were this effect to be significant, and 
also not to cancel as just suggested, we would expect a 
strong upwards evolution from z « to 2 which we do 
see, followed by a slow decline to higher redshift, which 
is certainly not reflected in the data. 

In short, the various methods and arguments all point 
clearly towards a significant evolution of the Lya space 
fraction with redshift. The main uncertainty affecting 
the precise absolute value of is probably due to 

statistical uncertainties in the LEs and to the simple ex- 
tinction correction applied to derive it, not possible Lya 
losses due to apertures. 
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5.2. The downwardly evolving escape fraction and the 
properties of the intergalactic medium 

Beyond a redshift of around 5.7, the measured value of 
/esc" begins to decline, although initially this decline is 
weak and the deviation from our best-fit relationship at 
redshift 6.5 is not sig nific ant. Adding the z = 6.5 point 
of lOuchi et all (|2010t l and lKashikawa et all (|2006h to our 
fit does not change the result. However the z = 7 point 
lies at just 8 %, and is around 2a below both the best- 
fit fe/ c a -z relationship (Figure [T]) and the predictions 
at this redshift based upon the f^J c a -EB-v relationship 
(Figure[3]). The z = 7.7 point formally takes the value of 
fe/ c a =0, and is presented with an extremely conservative 
error that is likely to be grossly overestimated (see § 12. 3[) . 
In comparison to z = 5.7, / C L S J" has declined by a factor 
of at least 2 by z — 7. We have so far attributed the 
increase in to an evolution in the dust content of 

galaxies, and it would be an extravagant departure from 
this evolutionary trend were ISM evolution to suddenly 
cause a sharp drop in at z > 6. Several other 

mechanisms are, however, naturally able to explain this 
break in the trend. 

5.2.1. Leaking ionizing radiation 

As discuss ed previous ly and by, f or exa mple , 

iBunker et all (|2010D and iBouwens et ail (|2009bf) . the 
LyC escape fraction at z ~ 8 must have been around 20- 
50 % in order to reionize the universe, depending upon 
the clumping factor of neutral hydrogen. Thus, as the 
galaxy population embedded in the reionization epoch 
evolves into the population observed at lower redshifts 
(~ 3), it must also transition through a phase of modest 
average f^/ c c (« 0.1-0.2). At these redshifts, measure- 
ments of f^/ c c are eme rging that do seem to be consis- 
tent with these values (|Iwata et al.l l2009: Vanzell aet al.l 
2010), which continue to evolve to lower values with de- 
creasing redshift (see iSiana et al.l I2010f ). Furthermore, 
since at z w 7 we are looking through the nearest edge of 
the reionization epoch into a partially neutra l Universe 
(as d etermined by quasar absorption studies, iFan et aT] 
2006), substantial LyC leakage must occur from the z ~ 7 
galaxies in order to complete reionization. 

If we set /^ c ~ at z = 5.7 and hold all the other 
properties of the galaxy population constant (i.e. no 
strong evolution of galaxy metallicity or IMF), this es- 
timate of fc/ c c ^ 30% at z m 7 would also reduce the 
nebular emission line spectrum to 70% of its value at 
z rj 6. This in itself would be sufficient to bring the pre- 
dicted value for f^/ c a within la of the measured value at 
z = 7. Thus, even in the redshift 7-8 domain we suggest 
that the the drop in the Lya LF could be attributed to 
the drainage of ionizing photons. 

5.2.2. Neutralizing the intergalactic medium 

As the IGM shifts from ionized to neutral, Lya pho- 
tons scatter in gas that immediately surrounds galaxies 
(|Miralda-Escuddll998t lHaiman fe Spaansl 1 1999ft . This is 
expected to manifes t as a drop in the observed Lya num- 
ber c ounts or LF (jRhoads fc Malhotral 120011 ; IHu et all 
2002), that tails much farther into the reionization epoch 
than absorption tests in quas ar spectra. Pre v iously 
iMalhotra fe Rhoadsl (|2006f ) and lKashikawa et all (|2006T ) 



have used the evolution of the Lya LF to look for such 
signatures of a neutral IGM transition but found conflict- 
ing results. However, the raw differential comparison of 
LFs only tests the ionized fraction if the evolution of the 
underlying galaxy population is u nderstood to an equal , 
or preferably better, level and iDiikstra et al.l (l2007f) 
showe d that the evolution reported by iKashikawa et all 
(2006) can, for example, be explained purely by the evo- 
lution of the dark-matter halo population. I n a similar 
vein to our own analysis, IStark et al.l (|2010bh have sug- 
gested the fraction of LBGs showing strong Lya emission 
to be a preferable signature of cosmic re-ionization to the 
evolution of the Lya LF. Further, from the lack of Lya 
line emission i n six o ut of seven z ~ 7 galaxy candidates 
I Font ana et al.1 (|2010f ) suggest that an increasingly neu- 
tral IGM is responsible for reversing the observed trend 
of the increasing fraction of strong emitters at redshift 
below ~ 6. 

By recasting the problem in terms of the Lya escape 
fraction, we remove the question of halo evolution from 
the problem - any halo mass function evolution is ac- 
counted for by the LBG population that is used to com- 
pute fcZ°- The drop in the Lya LF is also reflected by 
the fes C a -z diagram, quite securely by z — 7. If we hold 
the ISM properties and constant, we see that be- 

tween redshift 6 and 8 we need to suppress > 50 % of the 
Lya luminosity. However, what this means for the neu- 
tral gas fraction is much harder to infer since the fraction 
of photons that scatters in the IGM depe nds on the ex- 
act wavelength with which Lya is emitted (lHaimanll2002T; 
| Santosll2004HMalhotra fc. Rhoads 2004; Verhamme et al.l 
1200a IDiikstra fc Wvithell2010D . All we can say with re- 
liability is that the average effective optical depth seen 
by emitted Lya photons at z ~ 7 would be about 1. 

In summary the dip in the observed Lya escape frac- 
tion beyond a redshift of 6 seems to be real and, hold- 
ing all other galaxy properties constant, a loss of around 
50 % of Lya photons needs to be accounted for by 
z = 7 — 8. Current data can be equally well described by 
the galaxy population emitting this fraction of LyC pho- 
tons, and by Lya photons seeing an IGM optical depth 
(at the velocity of the emitting galaxy) of around 1. Ob- 
servational discrimination between the two scenarios will 
remain extremely challenging, but basically calls for fur- 
ther deep spectroscopic observations of the z — 7 — 8 
narrowband and dropout candidates, most likely requir- 
ing extremely large telescopes. 

5.3. Evolution of the dust content of galaxies 

So far we have been taking advantage of the fact that 
we have measurements of the dust extinction in our sam- 
ples of Ha and UV-selected galaxies. We have used this 
to infer the intrinsic star-formation rate density of the 
populations, and from there calculated / e ^° using Equa- 
tions [2] and [3l These Equations connect the quantity of 
/esc" wrtn Eb-v, via the ratio of the Lya- and UV- 
derived measurements of the star-formation rate density, 
Pi,. However, in §|4]we defined an alternative relation- 
ship between f^/ c a and Eb-v, based upon analyzing in- 
dividual galaxies, where we provide an empirical rela- 
tionship between these two quantities (Equation [4]) and 
relate them simply through coefficients. Thus we have 
four quantities (p* ya , p+ v , Eb~v, and /^°), that are 
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related by the various coefficients discussed in the previ- 
ous Sections. 

In all the previous Sections we have made use of the 
measured values of Eb-v but instead we could ignore 
this measurement, and invert the problem: use the ob- 
served Lya and uncorrected UV star-formation rate den- 
sities at a given redshift to estimate Eb-v , using Equa- 
tion Has a closure relation. Thus, substituting Equa- 
tion [4] into Equation [2] we can write: 



E 



1 



B-V 



0.4(A: A 



c Lya ) 



x logi 



•Obs 

r *,Lya 



•Obs 



(6) 



ya 



Out to z ps 6 we take the data compiled in Table [U 
and compute the observed SFRD from either Ha or the 
UV, depending on the redshift. We then use Equation [5] 
to estimate the sample-averaged Eb-v at each redshift, 
independently of the attenuation measurements them- 
selves. In short we ignore the fact that these Eb-v 
measurements have been made, and see if we can recre- 
ate them. We show the result as black data points in 
the upper panel of Figure HI with the actual measure- 
ments shown by the the small gray symbols. We then 
hypothesize that the dust content of the universe may de- 
crease exponentially, and adopting a function of the form 
Eb-v(z) = Cebv ■ e~xv(z/ zebv), we fit the coefficients 
Cebv = 0.386 and zebv — 3.42. Or, the e— folding red- 
shift scale for the Eb-v evolution is » 3.4. We show this 
relationship in Figure [4] with the thick red line. 

By performing this experiment we are throwing away 
observational information and the plot becomes some- 
what noisier, but nevertheless it resembles an inverted 
version of Figure [T] Fundamentally the plot shows a 
decrease in the Eb-v of galaxies as redshift increases, 
which is consistent with the measurements (Table [T] and 
gray points). This decrease in the dust content of galax- 
ies with redshift is already much-discussed in the litera- 
ture for LBGs at z ~ 2 —7, based upon a gradual bluening 
of the UV slopes (e.g. lHathi et alJl2008HBouwens et al.l 
l2009af) . At higher redshift, there is however a tendency 
for our new method to estimate higher Eb-v compared 
to the measurements obtained directly from the UV stel- 
lar continuum. In the upper panel of Figure |4] we also 
show the best-fitting relationships derived in the previ- 
ous Sections (black solid line) , where we take the redshift 
evolution of f^/ c a and use Equation|4]to convert to Eb-v 
using our best-fit coefficients - naturally this line almost 
perfectly reproduces the gray points. 

The red line (fit to these data) runs slightly flatter than 
the black one (combined fits from the previous sections) 
and suggests a slightly higher Eb-v, an d and therefore 
dust content, than me asured at the highest redshifts in 
IBouwens etafl (|2009a| ). At z ~ 3-6 it runs l ower, how- 
ever, than the measurements of lHathi et al.1 (|2008| ) who 
obtain slightly higher dust attenuations from LBG sam- 
ples. 

It is interesting to further investigate how the dust 
obscurations we derive c ompare with oth er estim ates. 
The SPH modell ing of iNagamine et al.l ([2008(1 and 
iDayal et al.l (|2009f ) already discussed in § 13.21 both pre- 
dict higher dust attenuations than measured in the z — 6 
dropout populations at Eb y=0.15. Detailed SEP mod- 
eling of z ~ 6-8 galaxies by Sch aerer fc de Barrosl (|2010l ) 
also suggest the presence of dust in some high- z LBGs. 
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Fig. 4. — Upper. The evolution of the dust content of galaxies 
with redshift. Black points show Eg_y derived from the raw ob- 
served (i.e. not dust corrected) star-formation rate densities in Lya 
and the UV (also Ha) using Equation [6] The gray points show 
the actual measured values which in general are well-reproduced 
by our new method. The black lines shows the predictions based 
upon the fcJi a -z and f^ a -Eg_y relationships derived in §[3] and 
§ |4] respectively. The red line shows the best fitting exponential 
function to these (black) data-points, and is shown to run slightly 
flatter, predicting more dust at higher redshifts. Lower: Eg_y 

measurements from the upper plot but translated into J^" us- 
ing our f^ a —EB—v relationship (Equation [4] and Figure[2}. The 
gray points show the same data as Figure [T] the black line shows 
the preferred fajj£*—Z power-law. The Figure demonstrates that we 
would have arrived at approximately the same conclusions, even if 
we had no measurements other than the Lya and UV luminosity 
functions. 

Here we estimate Eb-v~ 0.08 based upon the new 
methodolo gy. Similarly the se mi-analytical approach de- 
veloped in lBaugh et~atl (|2005t) find E B -v~ 0.1 at z > 3 
when examining the LBG population, which is certainly 
compatible with our estimates in the redshift 3-5 do- 
main. 

Since the empirical relationship derived between f^/ c a 
and Eb-v relates the two quantities directly, for com- 
pleteness we convert our _EB-v - redshift estimates to 
f^ a through Equation |4] This enables us to approx- 
imately re-create the main observational result of this 
article, Figure [TJ which we show in the lower panel of 
Figure |U Here we show the f^/ c a estimates derived in 
this Section with black shapes, with the original points 
from Figure [T] shown in gray. In short the difference be- 
tween the two sets of points is that in the gray ones the 
measured dust attenuation has been applied to the UV 
star-formation rate density in the computation of f^/ c a 
whereas in the black points, this quantity has been es- 
timated directly from the observed star-formation rate 
densities, using Equation [6] As with Figure Q] this shows 
feZ a increasing with redshift, but the actual estimates 
of the dust attenuation in the individual samples have 
not been used in the derivation of this Figure. The over- 
all trend of Figure Q] is maintained, although significant 
scatter has been added to the plot. It shows that even 
were no Eb-v measurements available, our main result 
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would have taken the same form and the overall trend 
would have been the same. 

6. SUMMARY 

We have compiled fifteen Lya luminosity functions 
from the literature between redshifts and 8 and in- 
tegrated them over homogeneous limits to obtain Lya 
luminosity densities. We have performed the same cal- 
culations with Ha emitting galaxies at z < 2.3, and with 
ultraviolet selected/dropout samples at z > 2.3, together 
with their extinctions due to dust. We subsequently used 
these dust-corrected luminosity densities to estimate the 
sample-averaged, volumetric Lya escape fraction (/^ a ) 
as a function of redshift. In summary, we show: 

• That / e L ^ a increases monotonically from the < 1 % 
level at z rts to around 40 % by redshift 6. Over 
this redshift range, the evolution can be well de- 
scribed by a power-law of the form fc/ c a (z) = 
C ■ (1 + z)t, for which we obtain coefficients of 
C* = (1.67±g;|2) x 1(T 3 ;£ = (2.57±g;l|). This rela- 
tionship predicts that f^/ c a should reach unity by 
a redshift of z = ll.l±g;|. 

• By combining samples of galaxies at redshift 2-3 
for which f^f c a and Eb-v have been computed, we 
derive a new empirical relationship between these 
quantities. This provides an effective attenuation 
law for Lya photons that includes not only the ef- 
fects of dust absorption, but also those of resonance 
scattering by neutral hydrogen. This new relation- 
ship takes the form f^ c a = C Lya ■ io-OA-E B _ v -k Lya _ 
where k\^ ya — 13.8 and C LyQ =0.445. 

• By combining our new f^/ c a -EB-v relationship 
with the measured dust content of (UV- and Ha- 
selected) samples in our study, we predict how f^/ c a 
should evolve with redshift, making no reference to 
Lya observations. Between redshift and 6.5, we 
find that this prediction is fully consistent with our 
measurements. Thus we are able to relate the up- 
wards redshift evolution of f^/ c a to the general de- 
crease in the dust content of the galaxy population. 
We discuss other effects that could mimic this trend 
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but ultimately find all of them to be implausible. 

• Beyond a redshift of 6 we see a drop in f^/ c a that 
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galaxy formation. 
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